Quantum phases of atomic Fermi gases with anisotropic spin-orbit coupling 
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We consider a general anisotropic spin-orbit coupling (SOC) and analyze the phase diagrams of both balanced 
and imbalanced Fermi gases for the entire BCS-Bose-Einstein condensate (BEC) evolution. In the first part, 
we use the self-consistent mean-field theory at zero temperature, and show that the topological structure of the 
ground-state phase diagrams is quite robust against the effects of anisotropy. In the second part, we go beyond 
the mean-field description, and investigate the effects of Gaussian fluctuations near the critical temperature. This 
allows us to derive the time-dependent Ginzburg-Landau theory, from which we extract the effective mass of 
the Cooper pairs and their critical condensation temperature in the molecular BEC limit. 
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I. INTRODUCTION 

The realization of BCS-Bose-Einstein condentate (BEC) 
evolution with two-component atomic Fermi gases has re- 
ceived tremendous attention in the last decade DjJ, [2|] . In these 
experiments, the tuning of attractive interactions permits the 
ground state of the system to evolve from a weak fermion at- 
traction BCS limit of loosely bound and largely overlapping 
Cooper pairs to a strong fermion attraction limit of tightly 
bound small bosonic molecules which undergo BEC. The 
main difference between the BCS-BEC evolution problem and 
the simple BCS theory is that the Cooper pairing is not al- 
lowed only for fermions with energies close to the Fermi en- 
ergy but is also allowed for all momenta. 

All of the early theoretical works were concentrated on the 
balanced Fermi gases, i.e. both components have the same 
number and mass, suggesting that the evolution is not a phase 
transition but a smooth crossover, and hence the name BCS- 
BEC crossover. This prediction was then found to be in very 
good agreement with all of the observations [1,2]. Motivated 
by the experimental success with balanced Fermi gases, many 
of the recent theoretical works were concentrated on imbal- 
anced (population, mass and/or dimension) Fermi gases, sug- 
gesting various phases and phase transitions between them. 
Some of these predictions were also confirmed by the fol- 
lowing experiments [3-6], opening the door for new studies 
on more complicated systems with the hope of finding exotic 
phases of matter. 

Arguably, one of the very promising new systems to investi- 
gate is the spin-orbit coupled (SOC) atomic Fermi gases J3,§]. 
This is mainly motivated by the very recent success in realiz- 
ing SOC atomic BEC HI O, and by a practical proposal for 
generating a SOC Fermi gas with 40 K atoms 111 111 . Some of the 
very recent results on these systems can be summarized as fol- 
lows. For the two-body problem, it has been found that a two- 
body bound state exists for some types of SOC, e.g. Rashba 
or Dresselhaus types, even on the BCS side (a s < 0) of a res- 
onance [12] with arbitrarily small a s —> 0~, where a s is the 
s-wave scattering length. For the many-body problem, it has 
been found for balanced Fermi gases that the SOC increases 
the single-particle density of states, which in return favors the 



Cooper pairing so significantly that increasing the SOC, while 
a s is held fixed, eventually induces a BCS-BEC evolution 
even for a weakly-interacting system when a s — > 0~ (v^vh . 
Similar to what happens in the usual BCS-BEC crossover pic- 
ture of a balanced Fermi gas with increasing l/a s , the BCS- 
BEC evolution with increasing SOC turns out to be a smooth 
crossover but not a phase transition. Therefore, the ground 
state of balanced Fermi gases with or without the SOC is a 
topologically trivial gapped superfluid. 

However, the BCS-BEC evolution with Rashba-type SOC 
is found to become a phase transition for population- 111814201 
and/or mass-imbalanced [21] Fermi gases. It has been found 
that the SOC counteracts the population imbalance, and that 
this competition tends to stabilize the uniform topologically 
nontrivial gapless superfluid phases against the phase separa- 
tion. In addition, topological phase transitions associated with 
the appearance of momentum space regions with zero quasi- 
particle/quasihole energies have been found, the signatures of 
which could be observed in the momentum distribution or the 
single-particle spectral function [ 19, 20]. 

The anisotropic (in momentum space) nature of the SOC is 
also found to stabilize exotic superfluid phases. For instance, 
in sharp contrast to the no-SOC case where only the gapless 
superfluid phase supports population imbalance, both the gap- 
less and gapped superfluid phases are found to support popula- 
tion imbalance in the presence of a Rashba-type SOC lll9ll20ll . 
Similarly, again in sharp contrast to the no-SOC case where 
only the gapped superfluid phase supports population balance, 
both the gapped and gapless superfluid phases are found to 
support population balance in mass-imbalanced SOC Fermi 
gases when the mass difference becomes large enough [21]. 

In this paper, we extend our recent works El El], and 
study the effects of anisotropic SOC on the phase diagrams 
of both balanced and imbalanced Fermi gases throughout the 
entire BCS-BEC evolution. We analyze both zero and finite 
temperature phase diagrams, and the paper is organized as 
follows. First, we review the noninteracting and interacting 
two-body problem in Sec. [TT1 and calculate the binding en- 
ergy of the two-body bound-state in vacuum. Second, we 
study the many -body problem in Sec. [Till where we derive the 
mean-field theory at zero temperature, and use it to analyze 
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the ground-state phase diagrams of imbalanced Fermi gases. 
Then, we investigate the Gaussian fluctuations in Sec. lIVI near 
the critical temperature, and calculate the effective mass of 
the Cooper pairs and their critical condensation temperature 
in the molecular BEC limit. Last, our conclusions are briefly 
summarized in Sec. [V] 



II. TWO-BODY PROBLEM 

Before presenting our new results for the many-body prob- 
lem, let us first introduce the model Hamiltonian and review 
some of the recent results for the noninteracting and interact- 
ing two-body problem. 

For the noninteracting SOC fermions, the two-body Hamil- 
tonian (in units of K = 1 = ks) can be written as 



e k,t + 

St. 



ek.j, - Sk z 



</'k 



(1) 



where ij>\. — [a\, ^, a\. J with a\. a (a k)(T ) creates (annihilates) 
a spin-cr fermion with momentum k = (k x , k y , k z ), ek.o- = 
k 2 /(2m) is the kinetic energy, and 5k x = OL x k x — ia y k y and 
= a zk z are the spin-orbit fields with {a x , a y , a z } > 0. 
The eigenvalues of this Hamiltonian matrix are 



£k, s = £k,+ + s\f (e k> - + Sk.) +\Skj. 



(2) 



where s = ± labels the helicity bands, and e^ s = (e k .-|- + 
sek,4,) /2 are the kinetic energy average and half of the kinetic 
energy difference of j" and I fermions. The corresponding 
eigenfunctions uj = (u* S ,U2 s ) are given by u\ tS /u2, s = 
S k± /{e K - + S kz -s V(ek,-+5 fc J 2 + |5 kx | 2 ] . Throughout 
this paper, we mainly consider four analytically tractable spin- 
orbit fields: (i) a x — a and a y — a z — corresponding to 
an equal mixture of Rashba- B22I1 and Dresselhaus-type 112 311 
SOC (ERD), (ii) a z — a and a x = a y — corresponding 
to a fully aligned SOC (FA), (iii) a x = a y = a and a z = 
corresponding to a purely Rashba- or Dresselhaus-type SOC 
(PRD), and (iv) a x = a y = a z = a corresponding to a fully 
spherical SOC (FS). We note that ERD- and FA-type SOC are 
essentially the same for balanced Fermi gases, but these cases 
differ substantially for imbalanced Fermi gases as we discuss 
below in Sec. Hill 

For the interacting SOC fermions, it has recently been 
shown that the strength of the attractive particle-particle in- 
teraction g > is related to the two-body binding energy 
e b < in vacuum via, 1/g = (1/2) J2k. s V( 2£ k,s + e tfc _e f>)> 
where e t h = ma 2 is the energy threshold for the two-body 
bound state. As usually done, the theoretical parameter g can 
be eliminated in favor of the experimentally relevant s-wave 
scattering length a s via the relation, 1/g = — mV / (Ana s ) + 
^ k l/(2ek,+), where V is the volume. 

The bound-state equation is analytically tractable for the 
four cases mentioned above. For the simplest ERD- and 
FA-type SOC, a two-body bound state with energy e& = 
— l/(ma 2 ) exists only when a s > 0, showing that the ERD- 
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FIG. 1 : (Color online) The binding energy ej < of the two-body 
bound state are shown as a function l/(maa s ) for the ERD-, FA-, 
PRD- and FS-type SOC. 



and FA-type SOC do not have any observable effect on the 
two-body problem. On the other hand, in the case of PRD- 
type SOC, a two-body bound state exists even for a s < lfl2ll . 
and its energy is determined by lfl4l[l9h 



maa s 



= J1 



ma 2 




(3) 



In the weak SOC limit, when ma 2 <C |e&|, this expression 
-l/(ma 2 ) — ma 2 + 2a /a s up to the leading or- 



gives e b 
der in a, 



and it recovers the usual result in the a — > limit. 



However, in the strong SOC limit, when ma 2 ^> |e&|, Eq. (0) 
gives eh fa —(Ama 2 / e 2 )e 2 ^ maa3 \ which is exponentially 
small as shown in Fig. Q] We also obtain ps —OAAma 2 
at unitarity JT3, EH]. Similarly, in the FS-type SOC, an even 
deeper two-body bound state exists for a s < 0, and its energy 
is given by 



2ma 2 



ma ± 



1 



Q 



Am 2 ai 



(4) 



< 



(a s > 0). This 
2 in the weak- and 



where + (— ) sign is valid for a 
expression reduces to e b ~ — n 
e b ps -l/(mOj) — 2ma 2 in the strong SOC limits, and to 
e b = —ma 2 at unitarity. 

Having shown that different types of SOC give rise to sig- 
nificant differences with regards to the existence of the two- 
body bound state and its binding energy, next we analyze the 
many-body problem where these differences also play an im- 
portant role. 



III. MANY-BODY PROBLEM 

Let us first consider noninteracting (g = or a s — > 0~) 
balanced (N-f = Aj_ = N/2 or /i-j- = /ij, = ji) Fermi gases at 
zero temperature. For this purpose, and throughout this paper, 
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we conveniently choose the energy (length) scale as the Fermi 
energy ep (Fermi momentum kp) of N /2 fermions such that 
N = k F V/(3w 2 ). 

It has been shown that increasing the SOC for a nonin- 
teracting Fermi gas leads to a change in the Fermi surface 
topology, when the number of fermions in the +-helicity band 
(N+) vanishes [13], where N s = J^k^it 1 ~ £ k,s)- This oc- 
curs when the chemical potential fi goes below the bottom 
of the energy band, i.e. when fj, = 0, or when a increases 
beyond a critical value (a c ). In some ways, this is similar 
to the usual BCS-BEC crossover problem, where the quasi- 
particle/quasihole excitation spectrum changes behavior as a 
function of increasing the scattering parameter l/(kpa s ) at 
fi = 0, i.e. its minimum is located at a finite (zero) momenta 
when fj, > (ji < 0). 

For the ERD- and FA-type SOC, we obtain \i = ep — 
ma 2 /2, and thus setting fi = gives a c = kp/m. For 
the PRD-type SOC, we obtain ■m^fm^i^ + 9ma 2 /2) = 
k F for > leading to [i m ep — 3ma 2 /2 in the 
weak SOC limit when ma 2 <C ep, and also obtain 
fx = 2k F /(3Trm 2 a) - ma 2 /2 for fj, < 0. Thus, set- 
ting /x = gives a c = [4/(3n)] 1/3 kp/m sa 0.75fcp/m. 
Lastly, for the FS-type SOC, we obtain 2m*Jm(n + 
2ma 1 )\j2[i + ma 2 — kp, leading to /i = — 3ma 2 /2 + 
m 3 a 4 /[4(2m 6 a 6 + k% + kpy/4m 6 a e + kp)} 1 / 3 + (m 6 a 6 + 



kp/2 + kpy/^a 6 + fc|/2) 1 / 3 /(2m). Thus, again setting 
fx = Ogives a c = (l/4) 1 / 3 fc F /m « 0.63k F /m. 

We emphasize that the aforementioned change in Fermi sur- 
face topology for the noninteracting balanced Fermi gases is 
not a quantum phase transition but a smooth crossover lfl3l 
[l9ll . Ad ding the interactions does not effect this picture 
much |[T3l4l6lfT9ll2lll . and the crossover in the Fermi surface 
topology again occurs at fi = 0. However, this is no longer 
the case for the interacting imbalanced Fermi gases 1118142111 . 
and next we show how this crossover picture becomes a quan- 
tum phase transition in the presence of a population imbalance 
within the self-consistent mean-field theory. 



A. Mean-Field Theory 

In the absence of a SOC and at low temperatures, it is well- 
established that the mean-field theory is sufficient to describe 
the Fermi gases both in the BCS and the BEC limits, and that 
this theory also captures qualitatively the correct physics in 
the entire BCS-BEC evolution 1!|,|2|]. Hoping that the mean- 
field formalism remains sufficient in the presence of a SOC, 
here we analyze the resultant ground-state phase diagrams. 

For this purpose, we use the mean-field Hamiltonian 



H 



■£4 





A* 



Sk ± 

-A* 





-A 

-£k,t + Sk z 

Sic 



A 





\ 



0\, 



(5) 



"£k4 — / 



where tp^ 
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j, k u k 1 , u,_k,f, 0— k,4.] denotes the fermionic 
operators collectively, £k,o- = £ k,er — AV, ar, d A = 
5(dk,t a -k4) is the mean-field order parameter where (• • ■ ) 
is the thermal average. The mean-field thermodynamic poten- 
tial can be written as ]8|] 



k.A 



1 + A k „ 



(6) 



where T is the temperature, A = {1,2,3,4} labels the 
quasiparticle/quasihole excitation energies E^.x, A"k,A = 
tanh[£; kiA /(2T)], and £ k , s = e k , s - fi s where /j, s = + 
s/zjj/2. Here, the quasiparticle/quasihole excitation ener- 
gies i?k,A are determined by the eigenvalues of the Hamil- 
tonian matrix given in Eq. (|5). Following the usual pro- 
cedure, i.e. dCl/d\A\ — for the order parameter and 
iV-f + sN± = —dVL/dfi s for the number equations, we obtain 



the self-consistency equations 
2|A| 1 dEk.x 
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Nf ± Ni = 



4^ 9IAI 

k,A 1 1 



(*k,A - 1) 



k.A 



1 ± 1 dE, 



k.A 



8fl± 



(*k,A - 1) 



(7) 
(8) 



These equations are the generalization of the mean-field order 
parameter and number equations to the case of an anisotropic 
SOC, and they are consistent with the known results in the 
appropriate limits El El IH (see Sec. HUB). 

As usual, we checked the stability of the mean-field solu- 
tions for the uniform superfluid phase using the curvature cri- 
terion SIM], which says that the curvature of D, with respect 
to |A|, i.e. 



d 2 n 

<9|Ap 



k,A L 



1 dE KX d 2 Eu X 



|A| 0|A| d\A\ 2 



(X K x - 1) 



J_ ( dEk >* 

~2T V a|A| 



(9) 
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needs to be positive, where Ik. a = sech 2 [_Ek,A/(2T)]. When 
the curvature d 2 il/d\ A| 2 is negative, the uniform mean-field 
solution does not correspond to a minimum of ft, and a 
nonuniform superfiuid phase, e.g. a phase separation, is fa- 
vored. It is known that the curvature criterion correctly dis- 
cards the unstable solutions, but metastable solutions may still 
survive. This may cause minor quantitative changes in the first 
order phase transition boundaries l20ll . 

B. Quasiparticle/Quasihole Excitations 

While the eigenvalues of the Hamiltonian matrix given in 
Eq. (0 do not acquire a simple analytic form for a gen- 
eral SOC, next we discuss three limits where the quasipar- 
ticle/quasihole excitation energies simplify considerably, al- 
lowing for further analytical investigation. 

1. Balanced Fermi gases with Sk ± 7^ and Sk z 7^ 

First of all for balanced (£k t — 6c 4.) Fermi gases, we 
obtain CMllIlSi 




where si = S2 = + and S3 = S4 = — , andpx = — (— 1) A , i.e. 
Pi — Pi — + an d P2 — Pi — — ■ Since all i?k.A have no zeros 
and are always gapped for all parameters in k space, we expect 
the BCS-BEC evolution to be a smooth crossover for balanced 
Fermi gases even in the presence of a SOC B13l4l6ill9ll . This 
is similar to what happens in the usual BCS-BEC crossover 
picture of a Fermi gas with no SOC, where the quasiparti- 
cle/quasihole excitation energies are also gapped. Therefore, 
the ground state of balanced Fermi gases with SOC is a topo- 
logically trivial superfiuid lfl9l - l2lll . 

In this case, the order parameter equation reduces to 
l/g = (1/2) ^ k Xk.s/(2£"k.s) and the number equa- 
tion reduces to N = (1/ 2) ^ J l - (£ k)S /.E k , s )X k , s ], 
where S k .+(-) = -Ek,i(2) Il3l4l7l l22tl . These equations 
are analytically tractable in the strong-coupling limit when 
fi < and 3> |A|, for which they are approximately 
given by l/g a (1/2) £ ks VPfck,. - m)] and N a 
(|A| 2 /2)X\ l/[2(£ k , s - A*) 2 ], respectively, where e k , s is 
given in Eq. Hj. Comparing the order parameter equation with 
the bound-state one immediately leads to [i = (e^ — ma 2 )/2 
for all four types of SOC, but the number equation leads to 
|A| 2 = 16y/e 3 F (2\n\ -ma 2 )/(37r) for the ERD- and FA-, 
|A| 2 = 2V2fci(2| fi\ — ma 2 )/(3nmy/m\iT\) for the PRD- 

, and |A| 2 = 2fc|(2|/x| - ma 2 ) 3 / 2 /(37rmV^|/i|) for the 
FS-type SOC. We note that |A| is independent of a only in 
the ERD- and FA-type SOC, which is consistent with the re- 
cent numerical findings lTi^[l6ll . This is not surprising since 
the SOC term can be eliminated by a momentum shift in the 
x and z directions, respectively, in the ERD- and FA-type 
SOC, which also leads to a shift in the chemical potential 
fi(a) = fi(0) — ma 2 /2. Since Eq. ( fTOt depends only on the 



total magnitude of the SOC, i.e. yJ\S kx \ 2 + S^, the ERD-, 

FA-, PRD- and FS-type SOC differ in their Jacobians of the 
k-space integrals. 

2. Sk x — > and Sk z 7^ 

When 5 ki -> and S kz ^ 0, i.e. the FA-type SOC, we 
obtain 

£ k ,A = s x y/fa,+ + P xS kz ) 2 + \A^ + PX &,-, (11) 

which again can be gapless in some k space re- 
gions. The zeros of _E k ,A can be found by impos- 
ing the condition E kM2) E K3{i) = ££_ - (£ k ,+ + 
P\Sk z ) 2 — |A| 2 = indicating that the zeros oc- 
cur at real (fcj. = \Jk 2 + k 2 , k z ) momenta such that 

k\ s = 2m/i + - k 2 z s + 2mS kz s + 2ms^ 2 _ - |A| 2 . 
Since k± s > 0, setting k± s = above leads to k z s — 

ma ± rJ 2m(/i + + ma 2 /2) + 2ms^J [i 2 _ — | A| 2 , provided 

that |A| < and |A| 2 < -(/x t + ma 2 /2)(n i + ma 2 /2). 
This analysis shows that the conditions |A| = \n—\ and 
|A| 2 = —(m + ma 2 /2)(ni + ma 2 /2) determine the phase 
boundaries between the SF, GSF(I) and GSF(II) regions (see 
Sec. IIII Q . and that these three phases meet at a tri-critical 
point determined by /i + = —ma 2 /2. 

In this case, the derivatives of the quasiparticle/quasihole 
energies are given by dE^x/9\A\ = s A |A|/ -B k , PA for the 
order parameter, dE K x/d h i+ = -s A (£ k:+ + P\S kz )/B Kpx 
for the average chemical potential, and dE^.x/ 'd\i- = —px 
for the half of the chemical potential difference. Here, 
B Kpx = ^J{i\,,+ +pxS kz ) 2 + |A| 2 . We again note for 
this case that the order parameter equation reduces to 
2|A|/ 3 = (l/4)£ k A (d£ k ,A/c>|A|)Xk,A, and the number 
equations reduce to Nf ± N± = (1/4) J2k aK 1 ±1 )/ 2 + 

{dE K x/dn ± )X K x\. 

Similar to the balanced case discussed in Sec. HUB 11 we 
note that FA-type SOC term again can be eliminated for the 
population-imbalanced Fermi gases, by a momentum shift in 
the z direction B24I1 . This also leads to a shift in the chemical 
potentials, i.e. n a (a) = /J, a (0) — ma 2 /2. Therefore, FA- 
type SOC does not have any observable effect on the phase 
diagrams, which is in sharp contrast with the ERD-type SOC 
as discussed next. 



3. S kz -> and S k± ^ 

On the other hand, when S k — > and 5 kj _ 7^ 0, it is 
straightforward to show that 

E K x = sxy/£ i+ + tl_ + |A|2 + |S k J2 + 2p x A k , (12) 

where A k = J$~Jtf~ + + |A| 2 ) + |^ k J 2 ^, + . Note in 
this case that i?k,A can be gapless at some points/lines in k 
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space. The zeros of Ek,A can be found by imposing the con- 



dition El m El X4) 

12 I o 12 



= - + |A| 2 - |S k ^ 
4\A\ 2 \S k± \ 2 = 0, indicating that both |S k J = and 



6c tfk \ + |A| =0 needs to be satisfied. Therefore, in the 

„ = 



PRD-type SOC, the zeros occur when ||19L|21[l k r = 
and at real k z momenta, k 2 s — 2m/i+ + 2ms</ju? 
provided that |A| < for /i+ > 0, and |A| 2 



|A| 2 , 



< tor ^ + > L), and |A| < — ^t^i 
for /i + < 0. Similarly, in the ERD-type SOC, the zeros occur 



when k x = and at real k p = ,/ k 2 



k 2 momenta, fc 2 s = 



2m/i + + 2ms^J fi_ — |A| 2 , provided with the same condi- 
tions as above. This analysis shows for both PRD- and ERD- 
type SOC that the conditions |A| = and |A| 2 = -^t^i 
determine the phase boundaries between the SF, GSF(I) and 
GSF(II) regions (see Sec. IHICt . such that these three phases 
meet at a tri-critical point determined by /i + = IU9I421I1 . In 
Fig. |2] typical excitation spectra E^a of the PRD-type SOC 
are shown for the SF, GSF(I) and GSF(II) phases, illustrat- 
ing the k-space topology of their gapped/gapless excitations 
as discussed in Sec. IIII Cl in great detail. 
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FIG. 2: (Color online) The excitation spectra Z?k,A of the PRD-type 
SOC are shown for (a) SF phase (a = 0.275k F /m and P = 0.25), 
(b) GSF(II) phase (a = 0.275fc F /m and P = 0.5), and (c) GSF(I) 
phase (q = 0.35kF/m and P — 0.5). These data correspond to the 
red cross marks in Fig.|4jb). 



In this case, the derivatives of the quasiparti- 
cle/quasihole energies are given by dE-^^x/ d\A\ — 
(1 + PA^ k _/A k )|A|/i? ki A for the order parameter, 

dE KX /dfi+ ' = -[i + px(£ t _ + |5 k J 2 )M k ]a,-/^ k ,A 

for the average chemical potential, and dEk.x/d/j,- = 
- [1 +P\ (£k,+ + 1 A 1 2 ) Mk] Ck, - /E K A for the half of the chem- 





no-SOC 


Balanced 


ERD 


PRD 


SF 


gapped 


gapped 


gapped 


gapped 


GSF(I) 


1 surface 


N/A 


1 ring 


2 points 


GSF(H) 


2 surfaces 


N/A 


2 rings 


4 points 



TABLE I: The topological classification of uniform superfluid phases 
are summarized, depending on the number of zero energy quasipar- 
ticle/quasihole excitation energy surfaces, rings or points in k-space. 



ical potential difference. Therefore, the order parameter equa- 
tion reduces to 2|A|/.g = (1/2) Ek, s ( as k, s /5| A|)X k , s , 
and the number equations reduce to Nf ± N± — 

(l/2)E k J(l ±1)/ 2 + {dE Ks /dv±)X Ks }, where 
£*,+(-) = E K1{2) ii. 



C. Ground-State Phase Diagrams 

There are three phases in the phase diagrams lfl9l - l2Tll . 
While the normal (N) phase is characterized by A = 0, the 
uniform superfluid and nonuniform superfluid, e.g. phase 
separation (PS), are characterized by <9 2 £!/<9|A| 2 > and 
d 2 il/d\A\ 2 < 0, respectively, when A / 0. Furthermore, 
in addition to the topologically trivial gapped superfluid (SF) 
phase, the gapless superfluid (GSF) phase can also be distin- 
guished by the momentum-space topology of its excitations. 
Depending on the number of zeros of E^ \ (zero energy re- 
gions in k space), there are two topologically distinct gapless 
phases. For the ERD-type SOC, we have GSF(I) where E^\ 
has one, and GSF(II) where i? k .A has two zero energy rings 
in k space. Similarly, for the PRD-type SOC, we have GSF(I) 
where £" k .A has two, and GSF(II) where _E kj A has four zero 
energy points in k space. The topological classification of uni- 
form superfluid phases are summarized in Table I. In Fig. [2] 
we show the excitation spectra E^ x of the PRD-type SOC for 
the SF phase inEJa), GSF(II) phase in^b), and GSF(I) phase 
in |2fc), illustrating k-space topology of their gapped/gapless 
excitations. These data correspond to the points indicated by 
the red cross marks in Fig.[4jb). 

In Fig. [3] we show the ground-state phase diagrams of (a) 
ERD- and (b) PRD-type SOC Fermi gases as a function of the 
population imbalance P = (Nf — N±)/N and the scattering 
parameter l/(kpa s ). Here, the dashed blue and dotted green 
lines correspond to |A| 2 = — /i^/i^ and |A| = \fi-\, respec- 
tively, and they mark the SF, GSF(I) and GSF(II) phase bound- 
aries. We note that the regions bounded by the solid black 
lines show instability toward a non-uniform superfluid phase 
(PS). Since our classification of distinct topological phases ap- 
plies only to the uniform superfluid region, the dashed and 
dotted lines shown within these regions are solely for illustra- 
tion purposes. 

We find that while the ERD-type SOC does not have 
any observable effect on the balanced Fermi gases (see 
Sec. IIII BTK it gives rise to a phase diagram with very similar 
topological structure as that of the PRD-type SOC. This is not 
very surprising since the main difference between the ERD- 
and PRD-type SOC is the Jacobians involved in the k-space 
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FIG. 3: (Color online) The ground-state phase diagrams of (a) ERD- 
and (b) PRD-type SOC Fermi gases are shown as a function of P = 
(JV-f — N_i)/N and l/(kpa a ). Here, we choose a = 0.15/cf/jti in 
(a) and a = (0.15/\/2)W m in (b), so that the magnitudes of the 
SOC are the same in both figures. The phase labels are described in 
the text (see Sec. IIII Cb . Note that the regions bounded by the solid 
black lines show instability toward a non-uniform superfluid phase 
(PS), and the dashed and dotted lines shown within these regions are 
solely for illustration purposes. 



FIG. 4: (Color online) The ground-state phase diagrams of (a) 
ERD- and (b) PRD-type SOC are shown as a function of P = 
(JV-f — Ni)/N and a at unitarity, i.e. when l/(fcfa s ) = 0. The 
phase labels are described in the text (see Sec. IIII Cl . Note that the 
regions bounded by the solid black lines show instability toward a 
non-uniform superfluid phase (PS), and the dashed and dotted lines 
shown within these regions are solely for illustration purposes. In (b) 
the red crosses mark a point from each phase for which the corre- 
sponding excitation spectra and momentum distributions are shown 
in Figs.[2]and[5] respectively. 



integrals. This must be contrasted with the FA-type SOC, 
which does not have any observable effect on the system even 
in the presence of a population imbalance (see Sec. IIII B 21 i. 
and therefore, its phase diagram is exactly the same as that 
of the usual population-imbalanced Fermi gases without the 
SOC. Since this problem is well-studied in the literature OlIII], 
we do not discuss it any further. 

Comparing Fig. [3] with the a — >• limit dO], it is clearly 
seen that both the ERD- and PRD-type SOC are counteract- 
ing the population imbalance. On one hand, this competition 
always tends to stabilize the GSF phase against the PS, and 
therefore, at any given P, the system eventually transitions 
to a stable SF or GSF by increasing a, no matter how small 
l/(fcfra s ) is. This is best seen in Fig. |4] where the phase di- 
agrams are shown as a function of P and a at unitarity when 
l/(kpa s ) = 0. On the other hand, we find that while both the 
ERD- and PRD-type SOC stabilize the GSF phase against the 
N phase for low P due to increased density of states lfT3l[l4ll . 
they destabilize the GSF phase against the N phase for high 
P. 



In sharp contrast to the a = case where only the gapless 
GSF phase can support population imbalance, one of the in- 
triguing effects of the ERD- and PRD-type SOC is that both 
the gapless GSF and gapped SF phases can support population 
imbalance when a ^ 0. This is possible due to the anisotropic 
nature of the SOC Fermi gases in k space. In fact, it has re- 
cently been shown for the mass-imbalanced Fermi gases that 
both the gapped SF and gapless GSF phases can support pop- 
ulation balance when a ^ 112 ill . This is again in sharp con- 
trast to the a = case where only the gapped SF phase can 
support population balance, and it is possible solely due to the 
anisotropic nature of the SOC Fermi gases in k space. 

The transition from GSF(II) to GSF(I) leads to a change 
in topology in the lowest quasiparticle band, similar to the 
Lifshitz transition in ordinary metals and nodal (non-s-wave) 
superfluids. However, the topological transition discussed 
here is unique, because it involves an s-wave superfluid, 
and could be potentially observed for the first time through 
the measurement of the momentum distributions of j" and J, 
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fermions The momentum distributions are readily 

available from Eq. (0, and we illustrate the typical T = 
distributions of the SF, GSF(I), and GSF(II) phases in Fig. \5\ 
The distributions are anisotropic in k space, which follows 
from the anisotropic structure of E^a- I n addition, while the 
distributions of SF phase do not show sharp features, those of 
GSF(I) and GSF(II) phases are exactly n k ,t = 1 and n k 4 = 
for k-space regions where k± = and k z < \k z \ < k z+ . 
Therefore, a major redistribution occurs for the minority com- 
ponent (n k jj at the topological phase transition boundaries. 
For instance, at the GFS(II) to GSF(I) transition boundary, the 
sharp peak that is present near the origin vanishes abruptly. 

Although this topological transition is quantum in its na- 
ture, signatures of it should still be observed at finite T, where 
the observables are smeared out due to thermal effects. While 
the primary signature of this topological transition is seen 
in the momentum distribution, single-particle spectral func- 
tion 12711 as well as some thermodynamic quantities such as 
the atomic compressibility would also show an anomaly at the 
transition boundary. 



«k.t 



«k4 




0.5 1 

k±/kj? 



FIG. 5: (Color online) The momentum distributions nk,<r of the 
PRD-type SOC are shown for (a) SF phase, (b) GSF(II) phase, and 
(c) GSF(I) phase. These data again correspond to the points indicated 
by red cross marks in Fig.[4{b). 

All of these results are obtained within the self-consistent 
mean-field approximation, which is known to be reliable for 
the entire BCS-BEC evolution only near T = 0. Since 
the fluctuations dominate the physics at finite T towards the 
molecular BEC limit j25ll . to emphasize further the effects 
of finite T, next we discuss the Gaussian fluctuations near 
the critical temperature T c , i.e. the temperature at which the 



mean-field order parameter A vanishes. 



IV. GAUSSIAN FLUCTUATIONS NEAR T c 

One way to go beyond the mean-field (or saddle-point) ap- 
proximation and study the Gaussian fluctuations is to use the 
imaginary-time functional integral approach l25l l26Tl . Us- 
ing this approach and expanding the order parameter field 
around A = 0, one obtains the fluctuation action Sq — 
0-/T) E q> „ £ _1 (q, i/„)|A(q, v n )\ 2 , where L-^q, v n ) is the 
inverse fluctuation propagator, A(q, u n ) is the fluctuation 
field, and v n = 1-kTti is the bosonic Matsubara frequency. 
It is a straightforward task to calculate the propagator 



g 8 ^ 

k,A ,A e 



A, 



k+5 



X, 



Ex 



k+f, A 



Ex 



k-#.A e 



1- 



Ck+a, + C k _q _ - (S kx+ s±.S* 



h.c.)/2 



(C? 



k+2,+ 



|5 k±+ ^l 2 )(^- f> - + l^x-^l 2 ). 

(13) 



where A D = {1,3} and A e = {2, 4} sums over odd and even A 
values, respectively, h.c. is the Hermitian conjugate, Ck, s = 
f k _ + sS k:c , and 



E K x = flA^/Kk,- +P\S kz )i + |5 ki | 2 +px&,+ (14) 

are the quasiparticle/quasihole excitation energies. The calcu- 
lation of the fourth order fluctuations is lengthy and straight- 
forward, but the results are not particularly illuminating. Next 
we use the fluctuation action to study the time-dependent 
Ginzburg-Landau functional near T c . 



A. Ginzburg-Landau theory near T c 

The Ginzburg-Landau theory is used to study the low- 
frequency and long-wavelength behavior of the order param- 
eter near T c . For this purpose, first we consider the static part 
of the propagator and expand L _1 (q, 0) in powers of qi, and 
then expand L _1 (0, v n ) — L _1 (0, 0) in powers of uj after the 
analytic continuation iv n — > ui + i0 + . 

This calculation leads to the time-dependent Ginzburg- 
Landau equation in k space rf25l, 26], i.e. L _1 (q, ui) = 
a {T) + Yli j CijQiQj/Z ~ dui. Here, the zeroth order coeffi- 
cient i _1 (0, 0) is given by 



1 1 



a m = ~-~ y ^ 

g 8 ^ S kA 

y k,A Q ,A 



1 - SA„PAe 



^k,A ~ ^k,A e 



^k,A e 



?k,- 



I5'k 



k,+ 



(15) 



The condition a(T c ) = is the Thouless criterion, and it leads 
to an equation for T c . We checked for all four types of SOC 
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that this criterion is in agreement with the order parameter 
equation after setting | A| = in the latter. The coefficient of 
the time-dependent term 



-X 



k,A ,A e 



1 ~ S\ PK 



k.A, 



<t2 



,) 



I 2 - si 



k,+ -r Pk_ 



2 )(Cf 



k,- 
Ek.X, 



-l^kj 2 ) 



k,A„ 



^k,A e ) 2 



CO 



(16) 



is a complex number. For balanced Fermi gases, while its 
imaginary part reflects the decay of Cooper pairs into the two- 
particle continuum for // > 0, the imaginary part vanishes for 
fi < and the behavior of the fluctuation field is propagat- 
ing, reflecting the presence of stable bound states (molecular 
bosons) lEME^ . 

The second-order coefficient d 2 L (q,0)/dqidqj evalu- 
ated at q = is not illuminating for imbalanced Fermi 
gases. However, for balanced Fermi gases with PRD-type 
SOC, where c xx = c yy = cj_, c zz = c z and c^j = 0, the 
coefficients are given by 



V 



' iJ 16m ^ I 2mT 2 £ , k 

k,s 



-kikj [ 1 



2TE k 
Els 



ma 

+ s l- 

ma kikj 
— 1£ 



Sij - 



ss 



Ek.* k 2 



k^ kj 

Sin 



(17) 



B. Molecular BEC limit 

For illustration purposes, here we consider only balanced 
Fermi gases. In the molecular BEC limit, when /1 < and 
\fi\ T c so that E^s/T — > 00, we may set X^ s — > 1 and 
^k,s — > 0. Therefore, in this limit, the time-dependent coeffi- 
cient simplifies to d = J2k s VC^-^k s) f° r a ^ types of SOC. 
The second-order coefficients simplify to c± = ^ k [1 + 
S ma/(2k ± )]/(WmE 2 s ) + £ k s 5 , SS '/[8k 2 ± (E Ks + E k , s ,)} 
along the (x, y) directions and to c z = ^ k l/(16m_E k ) 
along the z direction for the PRD-type SOC, and sim- 
ilarly to c = Ek.J 1 + 2sma/(3k)]/(16mEl iS ) + 
J2k s s> ss' /[6k 2 (E k . s + J5k, a ')] along all (a;, y, z) directions 
for the FS-type SOC. In the case of ERD- and FA-type SOC, 
we note that since (5^ — kikj/k 2 ) = for i = j = x and 
(Sij — kikj/k 2 ) = for i = j = z, respectively, the diag- 
onal coefficients all become equal in the molecular limit, i.e. 
c = Cu, and it is c = J2k,s V( 16m - E 'k, s )- 

These k-space sums are analytically tractable for all four 
types of SOC that we consider in this paper. For in- 
stance, for the ERD- or FA type SOC, we obtain d — 
2toco = my/mV/ (87ry2|/i| — ma 2 ). However, for the 
PRD-type SOC, we obtain d = my/2m\fj,\V/[8iv(2\n\ - 
ma 2 )] for the time-dependent, and c± — y/2mV(4\fj,\ — 
ma 2 )/[64TTy/\ir\(2\fi\ 



ma 2 )]-V2m/ [647r v /|7T|] ln[(2 \fi\- 
ma 2 ) / (2\fi\)] along the (x,y) directions and c z = 
m 2 i/2m|/i|l / /[47r(2|/^| — ma 2 )] along the z direction for the 
second-order coefficients. Similarly, for the FS-type SOC, 
we obtain d = my/rn\fi\V/[4Tr(2\fi\ — ma 2 ) 3 / 2 ] for the 
time-dependent and c = \/m(7\fj,\ — 3ma 2 )V/{24n(2\[i\ — 



v 2 W 2 l 



/(127ry / 2|/i|) for the second-order coeffi- 
cients. Next we extract the effective Gross-Pitaevskii param- 
eters using the asymptotic forms of the Ginzburg-Landau co- 
efficients given above. 



where E k + /_ 



+(-) 

a terms and the last J2k 



±i?k,i(2)- One needs to set the explicit 
, term to zero to extract c z from 
this expression. For balanced Fermi gases with FS-type SOC, 
where Cj, = Co and c^j = 0, one simply needs to replace fcj^ 
with k to extract cq. Similarly, in the case of ERD-type SOC, 
where c xx 7^ c yy — c zz and c^j = 0, one needs to replace 
k± with k x to extract c xx , and set the explicit a terms and the 
last ^ k , term to zero to extract c yy = c zz . Finally, in the 
case of FA-type SOC , where c xx = c yy 7^ c zz and c^j = 0, 
one needs to replace k± with k z to extract c zz , and set the 
explicit a terms and the last ^ k , term to zero to extract 

C-xx C yy - 

In general, the coefficients a(T), d and ca need to be calcu- 
lated numerically together with the order parameter and num- 
ber equations. However, it can be analytically shown that their 
asymptotic forms recover the usual Ginzburg-Landau equa- 
tion for BCS superfluids in weak coupling and the Gross- 
Pitaevskii equation for a weakly interacting dilute Bose gas 
in strong coupling 11251 12611 . Next we use the latter correspon- 
dence to extract the effective mass of the Cooper pairs (molec- 
ular bosons) and their critical condensation temperature in the 
molecular BEC limit. 



1. Effective Molecular Mass 

We recall that, since the Ginzburg-Landau theory derived 
above reduces to the Gross-Pitaevskii theory of a weakly- 
interacting molecular Bose gas after the rescaling ^(q, to) — 
v / dA(q, cj), the effective mass of the Cooper pairs (molecu- 
lar bosons) along the ith direction is simply given by ms,j = 
d/cii ll25ll26Tl . In the absence of a SOC, this gives = 2m 
for all (x,y,z) directions. 

Using the asymptotic forms of d and Co, and setting \n\ = 
(| £6 1 +ma 2 )/2, we find that the mass of the molecular bosons 
is itib,x — m B,y — nT>B,z = ttlb = 2m for the ERD- and FA- 
type SOC. However, using the asymptotic forms of d, cj_ and 
c z , and setting = (\ej,\ + ma 2 )/2, we find that the mass of 
the molecular bosons m,B x = tub 



2k| 



1^1 



■In 



2m 

iiib.i 2|efc| + 2mci 2 2|e&| + 2ma 2 \\eb\ + ma 2 

(18) 

and mB, z = 2m along the z direction for the PRD-type SOC. 
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FIG. 6: (Color online) The effective mass of the Cooper pairs are 
shown in the molecular BEC limit as a function l/(maa s ) for the 
PRD- and FS-type SOC. 



This expression is in complete agreement with the recent find- 
ings iThI [l5ll . showing that rriB.± decreases monotonically 
with increasing l/(maa s ) as plotted in Fig. [6] Equation (fl~8T l 

gives rriB = 4m when l/(maa s ) — » —00, ms = 2m 
when l/(mcra s ) —> 00, and tub ~ 2.40m at unitarity when 
l/(mcra s ) — > 0. On the other hand, we find that the mass of 
the molecular bosons rag x = rriB.y = rriB, z — mB is 



2m 
m B 



f6 



3/2 



2ma 2 



(19) 



for the FS-type SOC, which is also a monotonically decreas- 
ing function of l/(maa s ), as shown in Fig. [6] Equation ( fl9l i 
gives tub = 6m when l/(maa s ) — > —00, = 2m when 
l/(maa s ) — > 00, and m# = Z*/2m/ '(2y2 — 1) « 2.32m at 
unitarity. Having calculated the effective mass of the Cooper 
pairs, we are ready to calculate their critical BEC temperature. 



ith direction. In the absence of a SOC, this gives Tbec 
0.218e F JH. 



ft 0.16 




FIG. 7: (Color online) The critical BEC temperature of the Cooper 
pairs (molecular bosons) are shown as a function l/(maa s ) for the 
PRD- and FS-type SOC. 

Since we expect the Cooper pairs to progressively become 
weakly-repulsive in the molecular BEC limit, and using the 
effective mass of the Cooper pairs found in Sec. II VB 11 we 
find Tbec ~ 0.218e f for the ERD- and FA-type SOC, 
Tbec ~ 0.218{2m/m B ,i_) 2/3 e F for the PRD-, nndT BE c ~ 
0.218(2m/m B )e F for the FS-type SOC. In Fig. [7] we show 
Tbec as a function l/(maa s ). While Tbec is independent 
of a for the ERD- and FA-type SOC, we find that Tbec 
increases monotonically with increasing l/(maa s ) for the 
PRD- and FS-type SOC. For the PRD-type SOC, we find 
Tbec ~ 0.137eF when l/(maa s ) — > —00, Tbec ~ 
0.218e_F when l/(maa s ) — > 00, and Tbec ~ 0.193ei? at 
unitarity SB. Similarly, for the FS-type SOC, we find 
Tbec ~ 0.0726eF when \/{maa s ) -00, T B ec ~ 
0.218e_F when l/(mao s ) -t 00, and T B ec ~ 0.188e F at 
unitarity. 



2. Critical BEC Temperature 

The a dependence of the Cooper pair mass in the molecular 
BEC limit has a dramatic effect on the finite T phase diagram 
of the system. For atomic Bose gases, Tbec is determined 
from the number equation, N B = J2k V '(e ek - B/TBEC - 1), 
where N B is the number and e kiS = J2i={x, y ,z} k f I & m B.t) 
is the kinetic energy of atomic bosons, with rriB,i their ef- 
fective mass along the ith direction. This leads to Tbec = 
27r[nB/(^/n j ;ms i iC(3/2))] 2 / 3 in three dimensions, where 
71b = Nb/V is the density of bosons and £(x) is the Rie- 
mann zeta function with C(3/2) ~ 2.61. Setting = n/2, 
where n — N/V — k\/(Zir 2 ) is the total density of fermions, 
we obtain 

Tbec « 0.218 — ^ep (20) 

in three dimensions, where ep = k F /(2m) is the Fermi en- 
ergy and rriB,i is the mass of the molecular bosons along the 



V. CONCLUSIONS 

To conclude, we extended our recent works lfl9ll2lll . and in- 
vestigated the effects of an anisotropic SOC on the phase dia- 
grams of both balanced and imbalanced Fermi gases through- 
out the entire BCS-BEC evolution. We analyzed both zero 
and finite temperature phase diagrams, and our main results 
can be summarized as follows. 

In the first part, we derived the self-consistent mean-field 
theory at zero temperature, and used it to investigate the ef- 
fects of SOC on the ground-state phase diagrams. We showed 
that while both the ERD- and FA-type SOC do not have 
any observable effect on the balanced Fermi gases, only the 
FA-type SOC does not have any effect on the population- 
imbalanced gases. On the other hand, in the case of ERD- 
and PRD-type SOC, we found that the competition between 
the population imbalance and the SOC gives rise to very rich 
phase diagrams, involving normal, superfluid and phase sep- 
arated regions, and quantum phase transitions between the 
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topologically trivial gapped superfluid and the nontrivial gap- 
less superfluid phases. For instance, one of the intriguing 
effects of the SOC is that, in sharp contrast to the no-SOC 
case where only the gapless superfluid phase supports pop- 
ulation imbalance, both the gapless and gapped superfluid 
phases can support population imbalance in the presence of 
a SOC. We also showed that the topological structure of the 
ground-state phase diagrams is quite robust against the effects 
of anisotropy, i.e. they are very similar for ERD- and PRD- 
typeSOC. 

In the second part, we went beyond the mean-field descrip- 
tion, and investigated the effects of Gaussian fluctuations near 
the critical temperature. This allowed us to derive the time- 
dependent Ginzburg-Landau theory, from which we extracted 
the effective mass of the Cooper pairs and their critical con- 
densation temperature Tb ec m the molecular BEC limit. We 
showed that while the effective mass (Tbec) of the bosons 
does not depend on a for the ERD- and FA-type SOC, it de- 
creases (increases) monotonically as a function of increas- 
ing l/(maa s ) for the PRD- and FS-type SOC. We found 
Tbec « 0.14e F for the PRD- and T B ec « 0.073e F for the 



FS-type SOC in the weakly-interacting a s — > 0~ limit, and 
Tbec ~ 0.19e F for both types at unitarity when \a s \ — > oo. 
This shows that the presence of either a PRD- or FS-type SOC 
increases T c considerably especially in the BCS limit, which 
is in sharp contrast to ERD- or FA-type SOC where the SOC 
does not have any effect on T c throughout the BCS-BEC evo- 
lution. 
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